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Abstract 

The discovery of a light Higgs particle /z® (125 GeV) opens up new prospect for searching heavier Higgs boson(s) 
at the LHC Run-2, which will unambiguously point to new physics beyond the standard model (SM). We study the 
detection of a heavier neutral Higgs boson via di-Higgs production channel at the LHC(14TeV), 

WW*yY. This directly probes the Hhh cubic Higgs interaction, which exists in most extensions of the SM Higgs 
sector. For the decay products of final states WW, we include both pure leptonic mode WVT*—> ivtv and semi- 
leptonic mode WW*^ qq'tv. We analyze signals and backgrounds by performing fast detector simulation for the 
full process pp H ^ hh WW*-)'y (vivyy and pp H hh WW*yy —» ivqq'yy, over the mass range 
- 250 - 600 GeV. For generic two-Higgs-doublet models (2HDM), we present the discovery reach of the heavier 
Higgs boson at the LHC Run-2, and compare it with the current Higgs global fit of the 2HDM parameter space. 
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1. Introduction 

Most extensions of the standard model (SM) require an enlarged Higgs sector, containing more than one neutral 
Higgs states. After the LHC discovery of a light Higgs particle /z® (125 GeV) HE), a pressing task of the on-going 
LHC Run-2 is to search for additional new Higgs boson(s), which can unambiguously point to new physics beyond 
the SM. 

Such an enlarged Higgs sector 0 may contain additional Higgs doublet(s), or Higgs triplet(s), or Higgs sin- 
glet(s). For instance, the minimal supersymmetric SM (MSSM) 0 always requires two Higgs doublets and its next- 
to-minimal extension (NMSSM) Q further adds a Higgs singlet. The minimal gauge extensions with extra SU(2) or 
U(l) gauge group ||6l||7J will invoke an additional Higgs doublet or singlet. The minimal left-right symmetric mod¬ 
els IS] include an extra product group SU(2)^ ® U(l)£_^, and thus requires a Higgs bidoublet plus two Higgs triplets. 
For the demonstration in our present LHC study, we will consider generic two-Higgs-doublet models (2HDM) 19] 
under the SM gauge group. To evade constraints of flavor changing neutral current (FCNC), it is common to impose 
a discrete Z 2 symmetry on the 2HDM. For different model settings of Higgs Yukawa interactions, the 2HDMs are 
conventionally classified into type-I, type-II, lepton-specific, neutrino-specific, and flipped 2HDMs ||9l. The current 
study will focus on the conventional type-I and type-II 2HDMs. 

For the heavier Higgs state with mass above twice of the light Higgs boson h^, > 2M^ ^ 250 GeV, 

the di-Higgs decay channel H ^ hh is opened and becomes significant, in addition to the other SM-like major 
decay modes H VTVT, ZZ. Hence, the LHC can search for the di-Higgs production channel pp ^ H hh 
@|6l|IOllII]. ATLAS analyzed the decay channel hh bbyy at the LHC (8 TeV) run and found a 2.4cr excess 
at M{bbyy) ^ 300 GeV im. CMS performed similar searches for this channel and derived limits on the parameter 
space ca. An analysis of this channel at 14 TeV runs with high luminosity 1000 fb ' was done for 2HDM lEl. 
Another study considered the SM plus a heavy singlet scalar via // —> /z/z ^ bbWW* bbtviv channel for 14 TeV 
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runs with 3000fb“* luminosity ifTSll . We note that it is possible to increase the sensitivity of searches by studying 
and combining more decay channels of the di-Higgs bosons. 

In this work, we perform systematical study of production via a new decay channel of di-Higgs bosons, 
/7/7 —> H —> h/i —> . For the final state weak bosons, we will analyze both pure leptonic mode WW* tvlv 

and semi-leptonic mode WW* qq'iv. Since a SM-like Higgs boson /!°(125GeV) has decay branching fractions 
Br[/! —> bb, WW*, ZZ*] ^ (58%, 22.5%, 2.77%), we see that the di-Higgs decay mode hh —> WW*yy (with pure 
leptonic or semi-leptonic WW* decays) has the advantage of much cleaner backgrounds than hh —> bbyy, while 
Br[/! —> VkVk*] is only smaller than Br[/! —> bb] by about a factor of 2.6. Hence, we expect that the hh —» WW*yy 
mode should have comparable sensitivity to hh bbyy mode, and is more sensitive than hh —» bbWW* mode. 

This letter paper is organized as follows. In sectionl^ we present the production and decays of the heavier Higgs 
boson in 2HDM type-I and type-II. Then, in section]^ we systematically analyze the signals and backgrounds for 
the reaction pp ^ H ^ hh —> WW*yy, including both pure leptonic and semi-leptonic decay modes of the WW* 
final state. In section]^ we further analyze the LHC probe of the parameter space for 2HDM-I and 2HDM-II, and 
compare it with the current Higgs global fit. Finally, we conclude in section]^ 


2. Decays and Production of Heavier Higgs Boson //® in the 2HDM 


2.7. 2HDM Setup and Parameter Space 

For the present phenomenological study, we consider the 2HDM ||9] as the minimal extension of the SM Higgs 
sector. We set the Higgs potential to have CP conservation, and the two Higgs doublets Hj and have hypercharge 
T = -Hj , under the convention Q - + Y. It is desirable to assign a discrete Z 2 symmetry to the Higgs sector, 

under which the Higgs doublet IHIj (Hj) is ^^en (odd). With these, the Higgs potential can be written as 

V = H- + H^Hj) ^(H|Hi)2 h- ^(H^Hj)^ 

H-Tj IH 1 PIH 2 P -H T 4 IHIH 2 I 2 H- y [(HjH 2 )^ -I- (H^Hi)2] , (1) 

where the masses and couplings are real, and we have allowed a soft Z 2 breaking mass term of Mj 2 . The minimiza¬ 
tion of this Higgs potential gives the vacuum expectation values (VEVs), (EIj) = t'j)^ and (IHI2) = '^ 7 )^■ 

The two doublets jointly generate the electroweak symmetry breaking (EWSB) VEV v ^ 246 GeV, via the relation 
v-{v\ + v^^l^, where Vj = vcos/I and Thus, the parameter tan/I is determined by the Higgs VEV 

ratio, tany 6 = ^’ 2 /''! ■ Higgs doublets contain eight real components in total. 






(7= 1:2). 


( 2 ) 


Three imaginary components are absorbed by (VT*, Z“) gauge bosons, while the remaining hve components give 
rise to the two CP-even neutral states (hj, h^), one CP-odd neutral state A°, and two charged states 77*. The mass 
eigenstates {h, 77) of the neutral CP-even Higgs bosons are given by diagonalizing the mass terms in the Higgs 
potential Q. They are mixtures of the gauge eigenstates (/Zj, h ^), 

I h \ I cos a - sin O' \ / ^2 

\77 / \ sina cosaj\h^ 



where a is the mixing angle. Among the two neutral Higgs bosons, h is the SM-like Higgs boson with mass 
^ 125 GeV, as discovered at the LHC Run-1 HUll, and 77 is the heavier Higgs state. We will systematically 
study the LHC discovery potential of 77 state in the present work. The Higgs potential o contains 8 parameters in 
total, three masses and five couplings. Among these, we redefine 7 parameters as follows; the EWSB VEV v, the VEV 
ratio tanyS, the mixing angle a, and the mass-eigenvalues (M^,M^, M^, We may choose the 8 th parameter 

as the Higgs mass-parameter ■ Note that once we fix the mass spectrum of the 5 Higgs bosons as inputs, we are 
left with only 3 independent parameters (a, tanjS) and Mj 2 . The current LHC data favor the parameter space of the 
2HDM around the alignment limit ||9l, under which cos(j8 - a) - 0. This limit corresponds to the light Higgs state 
h to behave as the SM Higgs boson with mass 125 GeV. Eor practical analysis, we fix ^ 125 GeV by the LHC 
data and vary the heavier mass within the range of 250 - 600 GeV. We consider the Higgs states A and 77* to be 
relatively heavy, within the mass-range 47^,47^^ = 0.3 - 2TeV for simplicity. We will scan the parameter space and 
analyze the LHC production and decays of 77 in the next section. 
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Figure 1: Parameter space in - ( plane for 2HDM-I [plot-(a)] and 2HDM-II [plot-(b)], where the red dots present the viable points obeying 
the consistency requirement of the Higgs potential as explained in the text. 


The heavier neutral Higgs boson H has gauge couplings with Z‘^) and Yukawa couplings with quarks and 
leptons, which depend on the VEV ratio tan^S and mixing angle a. The gauge couplings of H with V {— W, Z) differ 
from the SM Higgs coupling by a scaling factor cos(J3-a), 

2M^ 

Gffyy = COS(/3-a)G^^yy, G^y - (4) 

The Yukawa interactions of H with fermions can be expressed as follows, 

f f - 

^y(H) = - Z (5) 

f=u,d,£ 

where the dimensionless coefficient differs between the Type-I and Type-II of 2HDM, as summarized in Table 1 


Table 1: Summary of the Yukawa couplings between the heavier Higgs boson and the SM fermions in 2HDM-I and 2HDM-II, where we 
have factorized out a common factor m^/v (corresponding to the SM Higgs Yukawa coupling). 


2HDM 

tu 

^H 

^H 

^H 

Type-I 

Type-II 

sino'/sinyS 

sino'/sin/? 

sino'/sinyS 

cosor/coSyS 

sino'/sinyS 

cosa/cosyS 


Inspecting the Higgs potential we derive the scalar coupling of trilinear vertex Hhh, 


^ Hhh 


cos(J3-a) 


cos(j3-a) 


[3MI -Ml- 2MI + [cos2O0-q') 


sin2(jS-a') 

tan2j8 


•Ml-V 


6Ml , ; 

- - 2MI 


^ sin2jS 


\ tan2y6 / sin2jS 


( 6 ) 


where in the second step we have used the relation = 2MIIsm2p. In the SM, the cubic Higgs coupling 

define a coupling ratio, ^ = G^jj^i^lG™!^ , which characterizes the relative strength of the Hhh 
coupling as compared to the h^ Higgs coupling of the SM. Under alignment limit cosfjS- a) —> 0, the trilinear scalar 
coupling (j^ takes the asymptotical form. 


Hhh 

/^sm 

^hhh 


(8M; 


j2/sin2/?- 


■Ml-2Ml) 


3Ml 


cos(j3-a) + 0{cos^(J3-a)). 


(7) 


In Fig.[^ we explore the parameter space of the Higgs potential 0 in the ^ plane. For ^ > 1, we expect 

that the decay branching fraction Br[// —> hh] and the production cross section cr[gg —> ^ hh] will be enhanced 
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Figure 2: Decay branching fractions of the heavier Higgs state for 2HDM-I [plot-(a)] and 2HDM-II [plot-(b)]. For illustration, we set tanyS = 1 

and = (500GeV, -(ISOGeV)^) for both plots. We also input cos(fi~a) = 0.4 for plot-(a) and cos(J3-a) = 0.1 for plot-(b). 


by the factor . In Fig.[^ the red points present the viable parameter space consistent with vacuum stability, unitarity 
and perturbativity bounds of the Higgs potential i). We also take into account the 3cr constraints from the current 
Higgs global fit (cf. Sec.|^. The electroweak precision data also constrain the parameter space of the 2HDM. It was 
found that in the 2HDM the charged Higgs mass satisfies, -600 GeV < < 100 GeV and > 250GeV 

ca, where is the mass of the heaviest neutral scalar. In the case with exact Z 2 {M -^2 - 0)> the potential could be 
valid up to the scale ~ lOTeV IfTTl . while for the present case of a softly broken the bound is much more relaxed, 
and the theory can be valid up to the Planck scale. For the analysis of Fig.[2 we have scanned the parameter space in 
the following range, tan/? e [1, 10], cos(J3-a) e [-0.6, 0.6], Mjj e [-200^, 200^] GeV^, e [200, 600] GeV, 
and 6 [300, 2000] GeV. In the following analysis, we will consider the same range of the 2HDM parameter 

space unless specified otherwise. 

2.2. Heavier Higgs Boson H®: Decays and Production 

Let us consider the decay modes of the heavier neutral Higgs boson H®. It is straightforward to infer the tree-level 
decay width for >2M„, 


r[H hh] 


327Tv^Mfj y 


( 8 ) 


For ^ 2M^ , we will include the off-shell decay H —> hh* with h* //, gg, yy , etc, where / denotes the light 
fermions except top quark. For the decay modes H —» VV, ff, we have Y{H —> yy]/F[// —> yy]sm = cos^(J3 - a) 
and F[// —> ff]/T[H —> //]sm = ■ (Here, the subscript “sm” denotes the “standard model” with a reference 

Higgs boson H which has the same mass as H in the 2HDM.) For the decay channel H ^ gg, we can express 
the partial width relative to the SM value, F[// —» gg]IT[H gg]sm = I ’ where = 

MjjKAmj:) and the function is the standard formula mM- The decay branching ratio of H ^ yy is 

practically negligible for Mh ^ 200 GeV. In Fig.[^ we present the decay branching fractions of the heavier Higgs 
boson H for both 2HDM-I [plot-(a)] and 2HDM-II [plot-(b)]. For illustration, we input tan/3 = 1 and - 

(SOOGeV, -(ISOGeV)^) for both plots. We also set cos(J3-a) — 0.4 forplot-(a) and cos(/3-a) = 0.1 forplot-(b). We 
see that for < 250 GeV, the dominant decay channels are H —> ZZ, WW, and for 250 GeV < < 350 GeV, 

the major decay channels include H —» ZZ, WW, hh since the H ^ hh channel opens up. For > 350GeV, the 
H tt channel is further opened, and will become dominant in 2HDM-II when cos(/3-a) takes values around the 
alignment limit as shown in Fig.|^b). But this situation can change when cos(yS-a) becomes larger and falls into the 
allowed region which separates from the alignment region (cf. Fig.|^in Sec.|^. 

From Eq. (j^ and Tabel[T] we see that the Yukawa coupling of the heavier Higgs boson H with tl has a scale 
factor = sina/sii\0 relative to the SM Higgs Yukawa coupling. The major LHC production channel is the gluon 
fusion process gg ^ H. Other production processes include the vector boson fusion pp Hqq', the vector boson 
associated production pp —> HV, and the top associated production gg —> Htt. The gluon fusion production cross 
section of H can be obtained from the corresponding SM cross section with a rescaling by // —> gg partial width. 


o-[gg^H] = (r[//^gg]/F[//^gg],Jo-[gg^//],^, 


(9) 
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Figure 3: Inclusive H production cross section via pp —» HX process at the LHC(14TeV), for 2HDM-I [plot-(a)] and 2HDM-II [plot-(b)] with 
tan/3 6 [1, 10]. All the red points satisfy the requirements of stability, perturbativity and unitarity, as well as the So* constraint by the current Higgs 
global fit. The cross section of inclusive H production pp HX contains four sub-channels from gg H, bh H, gb {gb) —» Hb (Hb), and 
gg {qq) —» Hbb. In plot-(c) and plot-(d), we present the sub-channel contributions to the inclusive cross section (T{pp HX^ for 2HDM-I and 
2HDM-II, respectively, where we set sample inputs, tanyS = 2 and cos(yS-Q') = -0.3 (-0.1) for 2HDM-I (2HDM-II). In plots (c)-(d), the red curve 
{gg HX contribution) and the black curve (summed total contribution) fully overlap because the gg HX channel dominates the inclusive 
cross section in the low tan/3 region. 


where we will include all NLO QCD corrections to the gluon fusion cross section as done in the SM case ESI. 
We note that for 2HDM-I, Table[2 shows that the H Yukawa couplings with top and bottom quarks have the same 
structure as in the SM, so the H production cross section <T[gg —> H] differs from the SM by a simple rescaling 
factor (sina/sin/?)^. For the 2HDM-II, we see that the H coupling to b quarks differs from that of t quarks by a factor 
of tauyS/tana, which can enhance the Z7-loop contribution to gg^H production for large tan/? region. Hence, the 
general relation (|^ should be used. The uncertainty of the gluon fusion cross section is about 10% over the mass-range 

= 250-600 GeV ED, which is roughly the total uncertainty of signal and background calculations. 

For the inclusive H production, we include the gluon fusion gg H, and Zj-related processes bh —> H, 
gb (gb) —> Hb (Hb), and gg (qq) —> Hbb. The production cross sections for these /j-related processes are derived by 
rescaling a factor of from their corresponding SM productions with the same Higgs mass. So we have the total 
inclusive cross section of pp HX for the 2HDM, 

o-[pp HX] = (F[// a-[pp(gg) -> H]^^ 

+ (^H? [o-[pp(bb) H]^^+ cr[pp(gb, gb) Hb, Hb\^+ cr[pp(gg, qq) Hbb],.J^ . (10) 

We present the inclusive H production rate for 2HDM Type-I and Type-II in Fig.[^a)-(b). Multiplying the produc¬ 
tion cross section with decay branching fraction Br(// ^ hh —> WW*jj), we compute the signal rate in the channej^ 
pp HX with H ^ hh ^ WW*jj. We summarize our results in Fig.for 2HDM-I and 2HDM-II, respectively. 
In Fig.|^a)-(b) and Fig.|2a)-(b), we have scanned the same 2HDM parameter space as in Fig.[2 The signal process 
is depicted by the left diagram of Fig.[^ From Fig.[^ we see that the cross section o-(pp —> HX) x Bt(H hh ^ 
WW*yy) can be as large as about 70 fb for 2HDM-I; while for 2HDM-II, this cross section can reach about 10 fb for 

< 340GeV. 


^Our analysis of the production rate of gg H ^ hh in the 2HDM is consistent with the recent study do). We thank Yun Jiang and Jeremy 
Bernon for providing data points of their calculation for numerical comparison. 
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Figure 4: LHC signal cross section cripp—>HX) xBt(H^ hh^WWyy) in the 2HDM with tan/3 £ [1, 10]. Plots (a) and (b) present the results 
for 2HDM-I and 2HDM-n, respectively. 


For comparison, we show the individual contributions of each sub-channel to the total inclusive cross section 
cr[pp HX] in Fig.|^c)-(d). For illustrations, we set sample parameter inputs, tanyS = 2 and cos(J3 -a) - -0.3 for 
2HDM-I, and tanyS = 2 and cos(j6 -a) - -0.1 for 2HDM-II. In plots (c)-(d), the red curve (gg —» H contribution) 
fully overlaps the black curve (summed total contribution). This is because the gluon fusion channel gg H 
dominates the inclusive production cross section for low tanyS region of the 2HDM. In general. Tab le[T] shows that 
for 2HDM-I, the H Yukawa couplings - sina/sin^S) are rather insensitive to tan/3. Hence, in 2HDM-I 

the gluon fusion actually dominates the H production over full range of tan j3 ^ 1, and the contributions of h-related 
sub-channels are always negligible. For 2HDM-II, the (up-type) H Yukawa coupling = sina/sinjS is the same as 
2HDM-I, and the down-type Yukawa coupling oc l/cos/3 - tan/I/sin/? is enhanced by a tan^S factor relative to 
. We find that for small tan /3 <3 , the gluon fusion channel still dominates the inclusive H production in 2HDM- 
II, and its cross section is larger than other h-related channels by a factor of (9(10-100) for ^ 300 GeV. The 
analysis of 2HDM-II in Sec.j^also concerns the small tan/3 region [cf. Fig.|^b)(d)]. Hence, in the following Sec.|3|^ 
we will focus our analysis on the Higgs production from gluon fusion channel, pp(gg) H ^ hh ^ 'WW*jy. 


3. Higgs Signal and Background Simulations 

In this section, we compute the Higgs signals and backgrounds at the LHC(14TeV). We perform systematical 
simulations by using MadGraphS package 11211 for the process, ppigg) H ^ hh ^ WW*yy, via gluon fusion 
channel. The parton-level Higgs production cross section cr{gg H) is derived from Eq. (|^, including NLO QCD 
corrections. We illustrate the signal Feynman diagram by the left plot of Fig.|^ For signal process, we generate 
the model hie using FeynRules Il22l . containing Hhh vertex and the effective ggH vertex. We compute signal and 
background events using MadGraphS/MadEvent ||2TI . Then, we apply Pythia Il23l to simulate hadronization of partons 
and adopt Delphes 11241 to perform detector simulations. 

Eor the hnal state WW decays, we will study both the pure leptonic mode WW —» Mv and the semi-leptonic 
mode WW —> qq'(v. The W decay branching fractions to ev and pv equal 10.8% and 10.6%, respectively, while 



diagram illustrates an irreducible background process gg qq'ivyy . 
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Table 2: Signal and background cross sections of pp WW*yy —» Mvyy and pp —> WW*yy —» qq'£vyy processes at the LHC (14TeV) after 
each set of cuts. The signal significance(Zo) is computed for the LHC(14TeV) runs with 300 fb“^ integrated luminosity. We input the heavier 
Higgs mass = 300 GeV, and set the sample signal cross section as (T{pp H hh —> WW^yy) = 5 fb. From the 3rd to 5th columns, we 
show the signals and backgrounds after imposing each set of cuts. The “Selection + Basic Cuts” are choosing according to Eqs. (ni-GD. In the 
pure leptonic mode, we impose the Final Cuts Mj{itw), M{£t), Mj{££vvyy), and AR{yy). In the semi-Ieptonic mode, we add 

the Final Cuts Py-(y) Mj{qq'£v), and A^(yy). 


pp —> {vlvyy 

Sum 

Selection-i-Basic Cuts 

Myy, Vt 

Final Cuts 

Signal (fb) 

0.525 

0.0251 

0.0214 

0.0161 

BG[^v^vyyH-^^yy] (fb) 

153.3 

0.937 

0.00225 

0.000215 

BG[ff/z] (fb) 

0.0071 

0.000493 

0.000419 

0.000076 

BG[Z/j] (fb) 

0.175 

0.0331 

0.00210 

0.000078 

BG[/!/!] (fb) 

0.00222 

0.000132 

0.000102 

0.000062 

BG[Total] (fb) 

153.48 

0.971 

0.00488 

0.00043 

Signihcance(Zo) 

0.734 

0.439 

3.70 

5.15 

pp ^ qq'{vyy 

Sum 

Selection-i-Basic Cuts 

Myy, Vj 

Final Cuts 

Signal (fb) 

2.2 

0.124 

0.0937 

0.0749 

^G[q^{vyy] (fb) 

31.59 

0.580 

0.0192 

0.00912 

BG[^vyy] (fb) 

143.3 

0.0642 

0.00349 

0.00182 

BG[W'/!] (fb) 

0.42 

0.00509 

0.00234 

0.00140 

BG[WW/i] (fb) 

0.0023 

0.000210 

0.000104 

0.000050 

BG[ff/z] (fb) 

0.0148 

0.00163 

0.000802 

0.000420 

BG[M] (fb) 

0.00462 

0.000291 

0.000160 

0.000106 

BG[f/i] (fb) 

0.0129 

0.000479 

0.000186 

0.000099 

BG[Total] (fb) 

175.35 

0.652 

0.0264 

0.0130 

Significance(Zo) 

2.87 

2.59 

7.29 

7.47 


that of W ^ TV is about 11.3% 11251 . The dijet branching ratio of W qq' equals 67.6% 11251 . Thus, the inclusion of 
semi-leptonic mode will be beneficial. Since t leptons can decay into e,/i, the detected hnal state e,ju will include 
those from the r decays. For M/, = 125 GeV, the branching fraction of h yy in the SM equals 2.3 x 10 m- 
In the following, we will first present the analyses for M„ - 300 GeV in Sec. 3.1 3.2 and then for heavier masses 


= 400,600GeV in Sec. 3.3 


3.1. Pure Leptonic Channel: hh y\/'W*yy^ tvtvyy 

For pure leptonic channel, we have hh —> 'W'W*yy Mvyy. Although this channel has an event rate about 
two orders of magnitude lower than that of hh bbyy mode, it has much cleaner background as compared to bbyy 
hnal state. After imposing simple cuts, we hnd that the backgrounds can be substantially reduced. We follow the 
ATLAS procedure for event selections. To discriminate the Higgs signal from backgrounds, we set up preliminary 
event selection by requiring two leptons (electron or muon) and at least two photons in the hnal state, 

rif = 2, tiy 2. ( 11 ) 

In the hrst step of event analysis, we need to prevent the potential double-counting, i.e., the reconstructed objects 
are required to have a minimal spatial separation ll26l . The two leading photons are always kept, but we impose the 
following criteria ll2^ : (i) electrons overlapping with one of those photons within a cone AR{e,y) < 0.4 are rejected; 
(ii) jets within A/?(jet, e) < 0.2 or A/?(jet, y) < 0.4 are rejected; (hi) muons within a cone of A/?(jU,jet) < 0.4 or 
AR{fi, y) < 0.4 are rejected. After this, we apply the basic cuts to take into account the detector conditions, which are 
imposed as follows. 


7’r(rX^r(?)>25GeV, PrW>15GeV, \q{y)\,\q{q)\Mf)\ < 2.5. (12) 

Next, we turn to the background analysis for pure leptonic mode. Besides the Mvyy and {{yy backgrounds, 
there are additional reducible backgrounds from Higgs bremstrahlung, vector boson fusion, and tth production. The 
cross section of the former two processes are fairly small and thus negligible for the present study. The tth associate 
production, with tt WWbb , can be important because the diphoton invariant-mass cut does not effectively discrimi¬ 
nate the signal process. But, this background can be suppressed by imposing h-veto Il27l . The production cross section 
for tth in the SM is cr(pp —» tth) - 0.6113 pb Il28ll . The latest h-veto efhciency of ATLAS is, e(h-veto) = 22% 
12^ . Thus, we estimate the cross section for this background process, 

cr(pp^tth^{v{vyy) - cr)/?/?—> ff/i) X Br[VT^^v]^ Br[/! ^yy] e(h-veto)^ ^ 7.28x 10^^ fb , (13) 
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where W ^ tv includes { - e,fi, t. We see that imposing the b-veto has largely suppressed the tih background. We 
note that the tth background is much smaller than the Ivtvyy background before kinematic cuts, while after all the 
kinematic cuts it could be non-negligible. So we will include both for the present background analysis. 





A(|)('Yyl 





AR(n) 



Mj.(llv\) (GeV) 



Figure 6: Signal and background distributions in the pure leptonic channel hh —> WW —> tvtv before imposing kinematical cuts. For comparison, 
we plot the signal distributions for = (300, 400, 600) GeV by (red, green, blue) curves. We present plot-(a) for Myy distribution, plot-(b) for 
Ifj distribution, plot-(c) for A^(yy) distribution, plot-(d) for AR(yy) distribution, plot-(e) for distribution, plot-(f) for AR(tt) distribution, 

plot-(g) for Mj{£tvv) distribution, and plot-(h) for Mj(t£vvyy) distribution, respectively. 
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Another potential background may arise from the Higgs pair production pp hh in the SM ll^lSTlI^ . Our 
signal process pp ^ H produces on-shell Higgs boson H with decays H hh, which has much larger rate as well 
as rather different kinematics from the non-resonant di-Higgs production in the SM. (Since our signal has on-shell H 
production, we find that its interference with the SM-type non-resonant hh production is negligible after kinematical 
cuts.) For instance, we can further suppress this SM di-Higgs contribution by imposing a cut on the transverse mass 
of di-Higgs bosons. 

We also consider a reducible background from the Zh associate production. The SM cross section of this process 
pp —> Zh at the LHC is o-{pp —» Zh) - 0.761 pb IfTFI . Hence, this background gives o-(pp Zh ^ i^jy) - 0.175 fb 
before any cuts. Because the Zh background must have the invariant mass M{i€) of final state di-leptons peaked at 

^ 91.2 GeV, we can efficiently kill this background by applying a narrow cut on which has little effect 

on the signal rate. In the present analysis, we choose, M{{€) e (M^-SF^, M^-i-SF^), where F^ - 2.5 GeV is 
the total width of Z boson. Other reducible backgrounds come from the fake events in which quark and/or gluon 
are misidentified as photons. These backgrounds include tvtvqy, Mvgy, Mvqq, Mvqg, and Mvgg. For our 
analysis, we adopt the fake rates used by ATLAS detector |[33l, 

- 3.6 X 10-^ - 3.6 x 10*^ (14) 

With such small fake rates, we find that these reducible backgrounds are negligible. 

In summary, with the above considerations of the SM backgrounds, we will compute the irreducible backgrounds 
with hnal state ivivyy , and the reducible backgrounds including the {€yy final state, the tth associate production, 
the Zh associate production, and the SM di-Higgs production. 

In Fig.|^ we present the distributions of relevant kinematical variables for the pure leptonic channel, including 
both signals and backgrounds. In this hgure, we show the signal distributions at the LHC(14TeV) with 300fb * 
integrated luminosity for = (300, 400, 600) GeV by (red, green, blue) curves as well as the backgrounds (black 
curves). Here we have input the sample cross section cr{pp H ^ hh ^ WW*yy) = (5, 3, l)fb for - 
(300, 400, 600) GeV, respectively. In the following, we will analyze how to effectively suppress the SM backgrounds 
by implementing proper kinematical cuts. 

From Fig.|^a)-(b), we hrst impose kinematical cuts on the diphotons invariant-mass M^y and the missing energy 
Ifj of final state neutrinos, 

120GeV <M^ < 130GeV, .^r>20GeV. (15) 

The missing energy cut can also sufficiently remove the €€yy background. 

Then, inspecting Fig.|^c)-(f), we apply the kinematical cuts on the azimuthal angle and opening angle A/? 
for the final state di-leptons and di-photons, respectively. 


A(p{U) < 2.0, ^R{{€) < 3.0, ^R{yy) < 3.8. 


(16) 


Here, from the distributions of Fig.|^c), we find that the /S.(p{yy) cut is not effective for Higgs mass - 300 GeV. 
So we do not implement this cut. 

For the transverse mass cut ll25l . we consider the transverse mass Mj for the ((vv system with two leptons and 
missing energy, which should be no larger than the Higgs mass ^125 GeV. All the final state leptons/neutrinos are 
nearly massless, so the transverse energy of each final state equals its transverse momentum Ej^ ^ \PT.i\, {i - 1,2,3), 
where i - 1,2 denote two leptons ^ and i - 3 denotes the system of two neutrinos. Thus, we have 

M\ = {Et,i + Et ,2 + - {Pt,i + Pt ,2 + Pt,z)^ - ^ 2£'j._i£'j.y(l-cos . (17) 

With this and inspecting Fig.j^g), we implement the transverse mass cut, 

Mjieevv) <\35G&W. (18) 

From Fig.j^h), we will further impose the transverse mass cut for the full final state ttvvyy , 


60GeV < MjiUvvyy) < 320GeV. 


(19) 


The kinematical cuts for the cases of - 400 GeV and 600 GeV will be discussed in Sec. 3.3 


We summarize the results in Table|^ for both signal and backgrounds. For demonstration, we first input the 
heavier Higgs mass = 300 GeV, and set the sample signal cross section cr(pp H^hh^ WW*yy) = 5 fb for 
the LHC (14TeV). In Table|^ we also show the significance of signal over backgrounds after each set of kinematical 
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Figure 7: Invariant-mass distribution of for semi-leptonic decay channel WW* qq'^v at the LHC(14TeV) with 300fb“^ integrated 
luminosity. Plot-(a) shows the mode with on-shell (olf-shell) decays qq' by green (blue) curve, for = 300GeV. The red curve 

corresponds to the realistic decays of WW* —> qq'lv. Plot-(b) presents the distribution for full signals of WW* —^ qq'iv by (red, green, blue) 
curves for = (300, 400, 600) GeV. In each plot, the black solid curve gives the full backgrounds. 


cuts at the LHC Run-2 with 300 fb ^ integrated luminosity. When the event number is small, we can use the median 
significance(Zo) (instead of 5/VB), as defined in following Il34ll . 


r 


Z„ = ^2 

(5+B)ln(__J-5 


As shown in Table|^ after applying all kinematical cuts, we estimate the signal significance(ZQ) = 5.15 . 


( 20 ) 


3.2. Semi-leptonic Channel: hh -^WW*yy^ qq'iv yy 


The analysis of semi-leptonic channel WW* qq'iv is similar to that of the pure leptonic mode WW* —> ivtv. 
But, there are nontrivial differences. One thing is that for each decay we need to specify which decay mode is from on- 
shell W (qq' or iv), since these two situations have different distributions. To illustrate this, we present the distribution 
of in Fig.j^a), where the green (blue) curve depicts the final state qq from on-shell (off-shell) W decays, and the 
red curve represents the actual distribution of from WW* —» qq'iv. Fig.lTla) shows that the distribution 
from on-shell W decays (green curve) has event rate peaked around - 70 - oOGeV, while the distribution 
from off-shell W decays (blue curve) is rather flat. 

Our first step here is also to remove the pileup events, similar to Sec. 3.1 Then, we select the final states by 
imposing the preliminary cuts 


nj^2, ^ 2 , rif-l. ( 21 ) 

For jets we choose the leading and subleading pair, while for photons we choose the diphoton pair whose is 
closet to = 125 GeV. Then, we choose the basic cuts to be the same as in Eq. ( fTS] ). 

Next, we turn to the background analysis. The most important background for this channel comes from the SM 
irreducible background, pp — > qq'ivyy, whose cross section is about cr[qq'ivyy] = 31.6 fb. Another significant 
reducible background is the SM process pp ivyy, which has a cross section cr[ivyy] = 143fb. But this will 
be mainly rejected by the jet-selections in Eq. ( |2T] l. Eor the tth background, we find that under b-veto its cross 
section is 0.0148 fb, as shown in Table]^ Single top associated Higgs production gives another background, cr[pp —» 
th(Ih) -I- A] = 79.4 fb 051 . where X represents single-jet or dijets in our simulation. We find that under b-veto this 
cross section of pp —» th(ih) -H A —> bivyy + X reduces to about 0.013 fb . We also include the non-resonant di- 
Higgs production in the SM, which has much smaller event rate and rather different kinematics. Other potential SM 
backgrounds may include the reducible backgrounds such as qqivgg with gg misidentified as yy. This is actually 
negligible due to the tiny g ^ y misidentification rate shown in Eq. ( [I4l i. 

Eor the kinematic cuts, we choose the M cut as in Eq. ( [T5| ). The invariant-mass should match the W 
mass. We depict the distribution in Eig.tTl Plot-(a) depicts the decay mode with on-shell (off-shell) decays 
W (W*) —> qq' by green (blue) curve, for = 300 GeV. The realistic decays of WW* —> qq'iv correspond to the 
red curve. In plot-(b), we present the distribution for full signals of WW* qq'iv by (red, green, blue) curves 
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for Mfj - (300, 400, 600) GeV. The black solid curve in each plot gives the full backgrounds. From Fig.|7j we choose 
the cut, 

M„„< 250GeV. (22) 




M-r (GeV) 






Figure 8: Signal and background distributions for semi-leptonic channel hh —» WW* —> qq'tvyy before imposing kinematical cuts. For com¬ 
parison, we plot the signal distributions for = (300, 400, 600) GeV by (red, green, blue) curves. We present the Myy distribution in plot-(a), 
the missing distribution in plot-(b), the Mj{qq'lv) distribution in plot-(c), the Pj{y) distribution of the leading photon in plot-(d), the A^(yy) 
distribution in plot-(e), and the Ai^iyy) distribution in plot-(f), respectively. 


We present the distributions for other kinematical observables in Fig.|^ where we have input the sample cross 
section (T(pp—>//—>/i/z—> WIT* 77 ) = (5, 3, l)fb for = (300, 400, 600)GeV. From Fig.[^a)-(b), we impose 
cuts on the diphoton invariant-mass Myy and the missing energy Ifj of final state neutrinos, 

120 GeV < Myy < 130 GeV, 10 GeV GeV. (23) 

We require Ifj > 10 GeV to suppress certain reducible backgrounds, as also adopted in the ATLAS analysis. For 
instance, consider the background qqjj+j with j mistagged as a lepton, where j denotes a gluon or quark jet. Since 
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Table 3: Signal and background cross sections of both pp —> 'W'Wyy —» Cvtvyy and pp —* VTW*yy —* qq'tvyy processes at the LHC(14TeV) 
after each set of cuts. The signal significance(Zo) is computed for the LHC (14TeV) runs with 300 fb“* integrated luminosity. We input the heavier 
Higgs mass =400GeV, and set the sample signal cross section cr(pp—>H—>hh-^WW'yy) = 3fb. From the 3rd to 5th columns, we present 
the signals and backgrounds after imposing each set of cuts. In the pure leptonic mode, we impose the Final Cuts Mj(Uvv), M{Ct), Mj(Uvvyy), 
A0(W). A(f>{yy), and Ai?('y')'). In the semi-leptonic mode, we add the Final Cuts F’j,(')'), Mj{qq'tv), A(^(yy), andAR(yy). 


pp —> Mvyy 

Sum 

Selection-i-Basic Cuts 

Myy, ll^T 

Final Cuts 

Signal (fb) 

0.315 

0.0165 

0.0147 

0.0107 

BG[(v(vyy+££yy] (fb) 

153.3 

0.937 

0.00394 

0.000169 

BG[fffi] (fb) 

0.0071 

0.000493 

0.000452 

0.000051 

BG[Zfi] (fb) 

0.175 

0.0331 

0.00247 

0.000065 

BG[/!/j] (fb) 

0.00222 

0.000132 

0.000116 

0.000074 

BG[Total] (fb) 

153.48 

0.971 

0.00698 

0.000359 

Significance(Zo) 

0.440 

0.289 

2.44 

4.05 

pp qq'tvyy 

^ total 

Selection-i-Basic Cuts 

Myy, ifj 

Final Cuts 

Signal (fb) 

1.32 

0.0891 

0.0671 

0.0533 

'&G[qq{vyy] (fb) 

31.59 

0.581 

0.0291 

0.00672 

BG[{vyy] (fb) 

143.3 

0.0642 

0.00454 

0.000891 

BG[Wh] (fb) 

0.42 

0.00509 

0.00335 

0.00139 

BG[WWh] (fb) 

0.0023 

0.000210 

0.000127 

0.000057 

BG[tth] (fb) 

0.0148 

0.00163 

0.00111 

0.000441 

BG[hh] (fb) 

0.00462 

0.000291 

0.000197 

0.000155 

BG[th] (fb) 

0.0129 

0.000479 

0.000247 

0.000104 

BG[Total] (fb) 

175.35 

0.653 

0.0386 

0.0098 

Significance(Zo) 

1.72 

1.87 

4.86 

6.22 


it contains no neutrino in the final state, we can eliminate it by imposing the missing energy Ijlj cut. This is more like 
a basic cut. For the transverse momentum distribution of the leading photon shown in Fig.[^d), we set the following 
cut, 

60GeV < Pj(y) < 150 GeV. (24) 

Then, we inspect the transverse mass distribution of qq'tv final state, which arises from the decay products of 
h —> WW*—^ qq'iv. From Fig.j^c), we impose the following cut, 

Mjiqqtv) < 200 GeV. (25) 

With Fig.[^e)-(f), we have also examined possible cuts on h.(p{yy) and AR{yy) distributions. We further impose, 

1 < ARiyy) < 3.8. (26) 

We summarize our results in Table|^ Here we present the signal and background cross sections after each set 
of cuts. We take an integrated luminosity of 300fb * for the LHC(14TeV), and derive the corresponding signal 
significance(Zo). Under all cuts, we estimate the final significance of the signal detection to be 7.47 in the semi- 
leptonic channel qq'tvyy, as shown in Table|^ 


3.3. Analyses of Heavier Higgs Boson with 400 GeV and 600 GeV Masses 


For signal and background analyses in Sec. 3.1 3.2 we have set the mass of heavier Higgs boson — 300 GeV 

for demonstration. In this subsection, we turn to the analyses for other sample inputs of Higgs mass, - 400 GeV 

and M„ - 600 GeV. We demonstrate how the analysis and results may vary as the Higgs mass increases. These are 


parallel to what we have done in Sec. 3.1 3.2 


For the heavier Higgs boson with mass = 400 GeV, from the distributions in Fig.|^ we choose the following 
kinematical cuts for the pure leptonic channel. 


120GeV < < 130GeV, A 0 ( 77 ) < 2.5, ^R{yy) < 2.5 , 

> 20 GeV, Mj{Uvv) < 135 GeV, 75GeV < Mjitvtvyy) < 420 GeV, (27) 

A0(^^) < 2.0, lARittf < 2.2, M{{€) \ M^+SV^) ■ 

Comparing with the previous case of = 300 GeV, we find that the distributions h.(p{{t), AR{tt), A(p(yy), 
and AR(yy) damp faster in the larger A(p and AR regions, as shown in Fig.[^ This is because the di-Higgs bosons 
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Table 4: Signal and background cross sections of both pp —» WWyy —> fvtvyy and pp —* WWyy —> qq'tvyy processes at the LHC(14TeV) 
after each set of cuts. The signal significance(Zo) is computed for the LHC (14TeV) with an integrated luminosity of 3 ab“*. We input the heavier 
Higgs mass = 600 GeV, and set the sample signal cross section <T{pp —>H—>hh^ WWyy) = 1 fb. From the 3rd to 5th columns, we present 
the signals and backgrounds after imposing each set of cuts. In the pure leptonic mode, we impose the Final Cuts Mj{CCvv), M{Ct), Mj(Uvvyy), 
A0(W), A(f>{yy), and Ai?('y')'). In the semi-leptonic mode, we add the Final Cuts Pj,(')'), Mjiqq'tv), Atpiyy), andAR( 7 ')'). 


pp —> tvtvyy 

Sum 

Selection-i-Basic Cuts 

Myy, ]}t 

Final Cuts 

Signal (fb) 

0.105 

0.00578 

0.00540 

0.00451 

^G[Mvyy+t£yy] (fb) 

153.3 

0.937 

0.00348 

0.000092 

BG[ff7!] (fb) 

0.0071 

0.000493 

0.000452 

0.000028 

BG[Z/!] (fb) 

0.175 

0.0331 

0.00138 

0.000029 

BG[/!/j] (fb) 

0.00222 

0.000132 

0.000117 

0.000070 

BG[Total] (fb) 

153.48 

0.971 

0.00543 

0.000219 

Significance(Zo) 

0.464 

0.321 

3.53 

7.76 

pp qq'tvyy 

^ total 

Selection-i-Basic Cuts 

Myy, ijlj 

Final Cuts 

Signal (fb) 

0.44 

0.0260 

0.0163 

0.0148 

BG[^^^vyy] (fb) 

31.59 

0.581 

0.00950 

0.00241 

P,G[tvyy] (fb) 

143.3 

0.0642 

0.00176 

0.000395 

BG[W/i] (fb) 

0.42 

0.00509 

0.00119 

0.000696 

^G[WWh] (fb) 

0.0023 

0.000210 

0.000035 

0.000035 

BG[ff7i] (fb) 

0.0148 

0.00163 

0.000402 

0.000237 

BG[/i/i] (fb) 

0.00462 

0.000291 

0.000120 

0.000087 

BG[f/i] (fb) 

0.0129 

0.000479 

0.000094 

0.000058 

BG[Total] (fb) 

175.35 

0.653 

0.0131 

0.00392 

Significance(Zo) 

1.82 

1.75 

6.70 

9.29 


are more boosted in the H ^ hh decays with heavier mass = 400 GeV . We present the cut efficiency for the 
case of Mj^ - 400GeV in Tablej^ where we set a sample signal cross section cr(ppHWW*yy) = 3 fb . In 
this case, we derive a signal significance(ZQ) = 4.05 after all the kinematical cuts. We also note from Fig.Qa)-(b) that 
in 2HDM-1 the cross section cr(/ 7/7 —> HX) x Br(// —> hh WW*yy) can reach up to 30 fb for = 400 GeV, while 
in 2HDM-II this cross section is below about 2fb at = 400 GeV . Hence, the significance for probing 2HDM-II 
with Mfj - 400 GeV will be rescaled accordingly, as we will do in Sec.|^ 

Then, we further analyze semi-leptonic channels for detecting the heavier Higgs boson H with mass - 
400 GeV . The corresponding signal and background distributions are presented in Fig.[^ Inspecting these distribu¬ 
tions, we choose the following kinematical cuts, 

120GeV < < 130GeV, < 250GeV, 

60 GeV < Pjiy) < 250 GeV, Mj{qq£v) < 250 GeV, > 10 GeV, (28) 

Aipiyy) < 2.3 , 0.75 < AR{yy) < 2.2. 

We summarize cut efficiency of the final state qqtvyy for = 400 GeV in Table|^ We derive a significance 
Zq = 6.22 after all kinematical cuts. 

Next, for the heavier Higgs H with mass = 600 GeV, the distributions of pure leptonic mode are shown in 
Fig.|^ From these, we set up the following kinematical cuts, 

120GeV < < 130GeV, ^T>25GeW, 

Mt{{(vv) < 135GeV, 75 GeV < MjiMvyy) < 620 GeV, (29) 

A4>{({) < 1.5 , AR({e) < 1.8 , M{ee) ^ (Mz-ST^, Mz+5Tz) , 

A(p(yy) < 1.8 , AR(yy) < 2.5 . 

The cut efficiency for = 600 GeV is summarized in Table|^ 

For the semi-leptonic final state qq'tvyy with = 600 GeV, we choose the kinematical cuts, 

120GeV <M^^< 130GeV, < 250GeV, 

Pj.(r) > 120GeV, Mjiqqtv) < 350GeW, > lOGeV, (30) 

A(p(yy) <1.6, AR(yy) <1.7. 
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With these, we summarize the cut efficiency of qq'tvyy final state for = 600 GeV in Tablej^ Since the typical 
production cross section with = 600 GeV becomes significantly smaller over the parameter space, we take a 
sample input cr(pp i/V)xBr(// ^ hh —> WW*yy) - Ifb, and consider an integrated luminosity of 3 ah ' at 
the LHC(14TeV). Hence, from Table|^ we can estimate the significance Zq = 7.76 and Zq = 9.29 for channels 
'W'W*yy Mvyy and WW*yy —> qq'ivyy, respectively. Besides, from Fig.|^a)-(b) we see that for = 
600 GeV, the cross section cripp —> //V)xBr(// ^ hh WW*yy) in 2HDM-I can reach up to 3 fb, while this cross 
section in 2HDM-II is below about 0.2 fb. Thus, the significance for probing the 2HDM-II with - 600 GeV will 
be rescaled accordingly. In the following Sec.[^ we will give a general analysis of the significance by scanning the 
parameter space of 2HDM-I and 2HDM-n without assuming a sample cross section. 

In the above analyses of Table[2]^ we have taken the sample cross sections, cr{pp H ^ hh ^ WW*yy') - 
(5, 3, l)fb, and an integrated luminosity Ji - (300, 300, 3000)fb“' for = (300, 400, 600)GeV. We have 
derived the significance of detecting H in each case. Thus, we may estimate the combined significance(Zo) by 
including both pure leptonic and semi-leptonic decay channels, 

ZQ(combined) = -^Z^iMvyy) + Z^iqq'ivyy) 

^ (9.06,7.41,12.1), for X = (300, 300, 3000) ft-'; (31a) 

^ (7.40, 6.05, 6.99), for £ = (200, 200, 1000)ft-'; (31b) 

which corresponds to = (300, 400, 600) GeV, respectively. 


4. Probing 2HDM Parameter Space at the LHC 

In this section, we study the probe of 2HDM parameter space by using the LHC Run-2 detection of the heavier 
Higgs state i/" via ppigg) -^H^hh^ WW*yy (Sec.[^, as well as the current global fit for the lighter Higgs boson 
/!''(125GeV) at the LHC Run-l. For the present analysis, we will convert the collider sensitivity (Sec.[^ into the 
constraints on the parameter space of 2HDM-I and 2HDM-IL As we showed in Fig.[^c)-(d) and explained in the last 
paragraph of Sec.|^ the inclusive Higgs production cross section a-{pp HX) is always dominated by the gluon 
fusion channel gg —> // in the small tanyS region, while other ^-related channels are negligible. (For 2HDM-I, this 
feature actually holds for full range of tany6 > 1.) Hence, for the present analysis, we will use Higgs production via 
gluon fusion ppigg) H ^ hh ^ WW*yy (Sec.|^ to probe the 2HDM parameter space. 

We combine the significance(Zo) from both pure leptonic channel WW*yy tvlvyy and semi-leptonic channel 
VTW*yy —» qq'tvyy at the LHC Run-2 with 300ft-' integrated luminosity. For this analysis, the relevant mass- 
parameters of the 2HDM are M^, . For demonstration, we will take the sample inputs, - 300,400GeV 

and (M^, = (500GeV, -(ISOGeV)^). With these, we have two remaining parameters in the 2HDM; the mix¬ 

ing angle a and the VEV ratio tan^S. In Fig.|^ we impose projected sensitivity of the LHC Run-2 by requiring 
significance(Zo) > 5 . From this, we derive the red contours in the parameter space of cos(jS-a')-tanyS plane, for 
2HDM-I [plots (a) and (c)] and for 2HDM-II [plots (b) and (d)]. The plots (a)-(b) correspond to = 300 GeV 
and plots (c)-(d) correspond to = 400 GeV. This means that the LHC Run-2 with an integrated luminosity 
-C = 300 ft"' can probe the red contour regions in each plot of Fig.with a significance(ZQ) > 5 . It gives a discovery 
of the heavier Higgs boson H (with 300 GeV or 400 GeV mass) in the red regions of the 2HDM parameter space. 

In Fig.|^ we further present the global fit for the lighter Higgs /!(125GeV) by using existing ATLAS and CMS 
Run-1 data, where the 2 cr and 3 cr contours of the allowed parameter space are shown by the green and yellow shaded 
regions, respectively. As we checked, our LHC global fit of the 2HDM is consistent with those in the literature li3^ . 
From this fit, we see that the parameter space favored by the current global fit is around the alignment limit of 2HDM 
with I cos(/l-a)| < 0.55 for 2HDM-I and | cos(yS-a)| < 0.15 for 2HDM-II. But, 2HDM-II still has an extra relatively 
narrow parameter region starting from tan /? > 2 . 

Fig.j^a) has input = 300 GeV for 2HDM-I. In this plot, the Zq > 5 region overlaps a large portion of 
the parameter space favored by the current LHC global fit. But, in Fig.j^b) for 2HDM-II, the situation is different 
because the overlap becomes smaller in the region coslyS-a) < 0, and gets enlarged for cos(y6-a) > 0. For the 
case of Mfj - 400 GeV in Fig.|^c), the probed parameter space of 2HDM-I has sizable reduction, especially for 
the region of cos(y6-a) > -0.05, in comparison with Fig.|^a) of - 300GeV . This is because the signal rate 
decreases as H becomes heavier [cf. Fig.|^a)]. On the other hand, for 2HDM-II, Fig.|^d) shows that the Zq > 5 
contours significantly shrink for = 400 GeV. This is because the signal rate for 2HDM-II drops more rapidly as 
Higgs mass rises to = 400 GeV in the small tanyS region [cf. Fig.^b)]. In this case, we see that the LHC Run-2 
with 300 ft-' integrated luminosity has rather weak sensitivities to the parameter space (shown by red contours), 
and the red contours no longer overlap with the favored region by the current LHC global fit (yellow and green 
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Figure 9: LHC probe of 2HDM parameter space in cos(^-a) - tan/3 plane. We impose projected sensitivity of the LHC Run-2 by requiring 
significance(Zo) > 5 for the process pp H hh WW*yy , with an integrated luminosity X = 300 fb“* at = 14 TeV. All red contours 
correspond to significance(Zo) = 5 with X = 300fb“^. Plots (a)-(b) [plots (c)-(d)] present the results for = 300GeV [M^ = 400GeV], 
while plots (a) and (c) [plots (b) and (d)] give the results for 2HDM-I [2HDM-II]. In plot-(d), the pink contours (Zq = 5) show a better probe with 
X = 3 ab“^ at the HL-LHC. The green (yellow) contours present the 2 <t (3(t) constraints from the Higgs global fit of 2HDM-I [(a) and (c)] and 
2HDM-II [(b) and (d)] at the LHC Run-l. In all plots, we have sample inputs (M^, M^ 2 ) - (SOOGeV, -(180GeV)^j, and the vertical dashed line 
denotes the alignment limit of the 2HDM. 


contours). We further analyze the probe from the upcoming High Luminosity LHC (HL-LHC) with 3 ah * integrated 
luminosity. We find that the HL-LHC can significantly extend the discovery reach of the parameter space of 2HDM-II, 
as demonstrated by the pink contour regions (Zp > 5) of Fig.|^d). 


5. Conclusion 

After the LHC discovery of a light Higgs boson /!‘’(125GeV) at Run-1, searching for new heavier Higgs state(s) 
has become a pressing task of the LHC Run-2. Such heavier Higgs state(s) exists in all extended Higgs sectors and 
can unambiguously point to new physics beyond the standard model (SM). 

In this work, we systematically studied the heavier Higgs boson production with the new decay channel, 
/7/7 ^ > /r/i —> WW* 77 , at the LHC Run-2. In section]^ we first analyzed the parameter space of the 2HDM 

type-1 and type-11, including the Hhh cubic Higgs coupling (Fig. [2. We computed the decay branching fractions and 
production cross section of the heavier Higgs boson H at the LHC Run-2 over mass range - 250 - 600 GeV, as 
shown in Fig.|2[Q Then, in section]^ we analyzed both pure leptonic mode WW* iviv and semi-leptonic mode 
WW* qq'tv. This channel has much cleaner backgrounds than the other process pp H ^ hh bbyy. We 
computed signal and background events using MadGraph5(MadEvent). We applied Pythia to simulate hadronization 
of partons and adopted Delphes for detector simulations. We followed the ATLAS procedure for event selections 
and built kinematical cuts to efficiently suppress the SM backgrounds. We analyzed various kinematical distributions 
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for pure leptonic and semi-leptonic decay channels in Fig.and Figs. [7j|^ for three sample inputs of Higgs mass 
— (300, 400, 600) GeV, respectively. In Table|2jQ we presented the signal and background rates of both channels 
under the kinematical cuts. In section]^ we combined the significance of pure leptonic and semi-leptonic channels, 
and analyzed the LHC Run-2 discovery reach of // as a probe of the parameter space in 2HDM-I and 2HDM- 
II (Fig.|^. For comparison, we further presented the current Higgs global fit of the LHC Run-1 data in the same 
plots. Finally, we note that it is hard to detect H with mass above 600 GeV at the LHC (14TeV) runs via di-Higgs 
production channel. We find it valuable to extend our present LHC study to the future high energy circular colliders 
pp(50-100TeV) 1371 , which are expected to further probe the heavier Higgs boson H with mass up to 6>(1 -5) TeV 
range via pp ^ H ^ hh production channel. 
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